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To develop an easy independent patient-specific quality assurance (QA) method using a benchmark plan 
for high-dose-rate intracavitary brachytherapy for cervix cancer, we conducted benchmark treatment plan- 
ning with various sizes and combinations of tandem-ovoid and tandem-cylinder applications with 'ideal' 
geometry outside the patient. Two-dimensional-based treatment planning was conducted based on the 
Manchester method. We predicted the total dwell time of individual treatment plans from the air kerma 
strength, total dwell time and prescription dose of the benchmark plan. In addition, we recorded the height 
(dh), width (dw) and thickness (dt) covered with 100% isodose line. These parameters were compared with 
169 and 29 clinical cases for tandem-ovoid or tandem-cylinder cases, respectively. With regard to tandem- 
ovoid cases, differences in total dwell time, dh, dt and dw between benchmark and individual plans were 
on average -0.2%±3.8%, -1.0 mm±2.6 mm, 0.8 mm±1.3 mm and -0.1 mm±1.5 mm, respectively. 
With regard to tandem-cylinder cases, differences in total dwell time, dh ftonl (the distance from tandem 
tip to tandem ring), dt and dw between benchmark and individual plans were on average -1.5% ±3.1%, 
—1.5 mm±4.9 mm, 0.1 mm± 1.0 mm and 0.2 mm±0.8 mm, respectively. Of two cases, more than 13% 
differences in total dwell time were observed between benchmark plans and the clinical cases, which turned 
out to be due to the use of the wrong source position setting. These results suggest that our method is easy 
and useful for independent verification of patient-specific treatment planning QA. 

Keywords: independent verification; treatment planning; Manchester method; benchmark plan; 
high-dose-rate intracavitary brachytherapy; uterine cervix 



INTRODUCTION 

Brachytherapy is an essential component of radiotherapy for 
the carcinoma of uterine cervix and is often combined with 
external beam radiation therapy (EBRT) for radical treat- 
ment. Several studies have suggested that control rates are 
significantly improved with EBRT and brachytherapy [1, 2]. 
High-dose-rate (HDR) remote afterloading intracavitary 
brachytherapy is widely used throughout Asia and Europe 
[3], and is becoming steadily more common in the USA [4]. 



The importance of independent verification of dosimetry 
prior to HDR brachytherapy treatment delivery has been 
recognized worldwide, and is specified in the guidelines 
of international regulatory agencies [5]. The Nuclear 
Regulation Commission (NRC) considers a 20% difference 
between the prescribed total dose and delivered dose to be 
a reportable medical event [5]. Thomadsen et al. identified 
44 medical events in HDR brachytherapy between 1980 
and 2001 in data from the NRC and International Atomic 
Energy Agency [6]. In fact, patients are often required to 
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wait during treatment planning with an applicator inserted 
by a radiation oncologist, during which time errors and 
miscommunications can easily occur. This situation clearly 
indicates that patient-specific quality assurance (QA), in- 
cluding independent verification of treatment planning 
and confirmation of applicator geometry, should be done 
quickly and easily. 

Many studies have reported independent verification 
methods for HDR brachy therapy treatment planning [7-16]. 
More recent reports have focused on the development of 
in-house software based on the AAPM TG 43 [17] formu- 
lism to calculate the dose at arbitrary points [13-15]. 
Although such software might be useful for the commis- 
sioning of treatment planning systems, human errors in in- 
dividual treatment planning in clinical practice will not be 
identified due to the use of the same coordinate system, 
digitized applicator paths and dose point coordinates as 
those in the treatment planning system. 

Although image-guided intracavitary brachytherapy has 
been enthusiastically investigated [18], treatment planning 
based on the Manchester method using two projection 
radiographs is still used [3]. One of the goals in intracavi- 
tary brachytherapy for carcinoma of the uterine cervix is to 
achieve the same level of consistency as the Manchester 
method. We have established the Osaka University Protocol 



based on the Manchester method with some modifications 
[19]. The goal of our institution is to achieve consistency 
with our protocol-based benchmark plans. 

Here, we propose a very quick, simple and easy patient- 
specific independent verification method for Manchester 
method-based treatment planning using benchmark plans to 
detect human errors and evaluate the quality of the applica- 
tor geometry in the patient. 



MATERIALS AND METHODS 
Creation of benchmark plans 

In this study, Fletcher-type (Fletcher- Williamson 
Asian-Pacific) tandem-ovoid and tandem-cylinder metal 
applicators (Nucletron International B.V., Veendaal, the 
Netherlands) were used. Various sizes and combinations of 
these applicators were constructed by one radiation oncolo- 
gist with the 'ideal' applicator geometry outside the patient 
(Fig. la and b) and then reviewed by a medical physicist. 
We constructed eight kinds of tandem-ovoid and six kinds 
of benchmark plans (Table 1). 

Figure lc shows the 'ideal' geometry of a tandem-ovoid 
applicator used in our institution. Namely, a flange on the 
tandem tube is used at the origin, which is the cervical os 





(c) 




Cervical os level 



Point A 




Cervical os level 



Fig. 1. Creation of benchmark plans, (a) Scheme for the construction of tandem-ovoid; (b) tandem-cylinder 
applications; (c) typical dose distribution with the tandem-ovoid application. 
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Table 1. Applicator settings used in benchmark plans 



Tandem-ovoid 



Tandem 
length (cm) 



Ovoid size 



Tandem-cylinder 



Cylinder diameter (cm) 



4 


S 


ss 


2 


2.5 


3 


5 


S 


ss 


N/A 


N/A 


N/A 


6 


s 


ss 


2 


2.5 


3 


7 


s 


ss 


NA 


NA 


NA 



The ovoid diameter of S size is 2.0 cm. The size of SS ovoid 
is half-cut-size of S ovoid. NA, not applicable. 

and the tip of the ovoid is aligned with the origin (a flange 
on the tandem). Using these ideal applicator settings, treat- 
ment planning was performed with PLATO (Nucletron 
International B.V). Source dwell time was manually opti- 
mized based on the Manchester method as demonstrated by 
Tod and Meredith with minor modification [19-21]. Air 
kerma strength and total dwell time were then recorded. All 
benchmark plans were constructed by a medical physicist 
and reviewed by another medical physicist. 



D(r, 6), x S KJ , x dwell time b 

Dwell ttme t = —, ; 3 

D{r, 0) h x S KJ 1 ' 

where D(r, 6) t and D(r, 0) h represent the prescription dose 
of each treatment and the benchmark plan, respectively. 



Comparison of dose shape with that of the 
benchmark plans 

ICRU report 38 [23] recommends reporting the reference 
volume as well as total reference air kerma strength and 
absorbed dose at reference points. The reference volume is the 
volume encompassed by the reference isodose surface, which 
is represented by the major dimensions of the following: 

(i) Height (dh), which is the maximum dimension 
along the tandem source measured on an 'oblique' 
sagittal plane; 

(ii) Thickness (dt), which is the maximum dimension 
perpendicular to the tandem sources measured on 
a transverse plane; 

(iii) Width (dw), which is the maximum dimension 
perpendicular to the tandem sources measured 
on a transverse plane. 



Prediction of total dwell time for individual patient 
treatments 

Dose was calculated with the following formula introduced 
in AAPM TG43-U1 [22]. 

D(r, 9)=S K xAx ^f' x g L (r) x F(r, 6) 

x dwell time 
= St x dwell time x A 

where S K , A, G L , g L and F represent air kerma strength, 
dose rate constant, geometric function, radial dose function 
and anisotropy function, respectively. Here, we defined A 
as the product of A, G L , g L (r) and F (r, 0). 

The dose at the reference point in the benchmark plan is 
therefore calculated by the following: 

D(r, 6) b = A h x S K _b x dwell time;, (1) 

Similarly, the dose at point A in the individual plan is cal- 
culated by the following: 

D(r, 6), = A h x S K _, x dwell time, (2) 

If the 'ideal' tandem-ovoid geometry is achieved in the 
patient without any planning errors or misuse of the appli- 
cator, A b is nearly equal to A t . Dwell time in the individual 
plan can therefore be predicted by the following formula 
from (1) and (2): 



In addition to the above parameters, we defined the 
dimensions of dh front and dh ext , which represent the dis- 
tance from th 100% isodose line of the tip side of the 
tandem to the origin and that from the 100% isodose line 
of the connector side of the tandem to the origin, respect- 
ively (Fig. 2c). 

Figure 2 shows the definitions of these parameters. For 
the tandem-cylinder, we recorded the additional parameters 
of dh front and dh ext , which represent the maximum dimen- 
sion of the 100% isodose line of the tip side of the tandem 
to the tandem flange, and that of the connector side of the 
tandem to the tandem flange (Fig. 2c). These values were 
measured for individual treatment plans and compared with 
those of the benchmark plans. 

Analysis of clinical cases 

We retrospectively analyzed 168 and 29 clinical cases from 
2009 through 2010 with a tandem-ovoid and tandem- 
cylinder, respectively. The difference in total dwell time 
between a benchmark plan and an individual treatment plan 
was calculated using the following formula: 

Relativedifference(%) = dividual - T benchmark ^ m 

"^benchmark 

where T benchmark and T individua i represent the total dwell 
time of the benchmark and individual plans, respectively. 
Differences in dose distribution shapes, including dh or 
dhf ront , an d dh ext , dt, and dw between the benchmark and 
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(a) (b) 




Fig. 2. Definitions of dh, dhfront, dhext, dt and dw used in this study based on ICRU report 38 [23]. These 
parameters were defined in the plane of the cervical os. (a) Transverse plane of tandem-ovoid case; (b) 
sagittal plane of tandem-ovoid case; (c) sagittal plane of tandem-cylinder case. Red line shows the isodose 
line of 100% of the prescription dose. 



individual plans were calculated using the following 
formula: 

Relativedifference = shape be nchmark - shape indiv i dua i 

Where shape benchmark and shape individua i represent the dose 
distribution shapes of the benchmark and individual treat- 
ment plans, respectively. 

Correlations between the differences in total dwell time 
and those in dh, dt or dw among the benchmark and indi- 
vidual treatment plans were evaluated by Spearman's rank 
correlation coefficient using Dr. SPSS II software (1MB, 
New York, USA). 

Tolerance levels for total dwell time, dh or dh front , and 
dh ext , dt, and dw were calculated by the following formula, 
which was first proposed by Venselaar et al. [24] 

Tolerance level = mean deviation ± 1 .96 SD 
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Fig. 3. Histogram of % differences in total dwell time between 
the benchmark and individual treatment plans. 



Differences averaged -0.2% ±3.8% (range, -13.3-9.6%), 
and exceeded 5% in 23 of 169 clinical cases. 



RESULTS 

Tandem-ovoid cases 

Prediction of total dwell time for individual 
treatment plans 

Figure 3 shows a histogram of differences in total dwell 
time between benchmark and individual treatment plans. 



Comparison of the dose distribution shapes of 
individual treatment plans with those of the 
benchmark plans 

Figure 4 shows a histogram of differences in dose distribu- 
tion shapes between the benchmark and individual treat- 
ment plans in tandem-ovoid applications. The differences in 
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Fig. 4. Histogram of differences in (a) dh, (b) dt and (c) dw between the benchmark and individual treatment 
plans in mm. 



dh, dt and dw between the benchmark and individual plans 
averaged -1.0 mm ± 2.6 mm (range, -8.6 mm to +6.5 mm), 
0.8 mm± 1.3 mm (range, -1.4 mm to +5.2 mm) and 
-0.1 mm ±1.5 mm (range, -5.1 mm to +4.0 mm), 
respectively. 

Regarding dh, 9 of 169 cases showed a difference 
between the benchmark and individual plans of greater than 
5 mm. For dt and dw, in contrast, only one case showed a 
difference of more than 5 mm. 

Subset analysis of cases with large deviations 
between the benchmark and individual treatment 
plans 

We verified that all tandem-ovoid treatment plans were ap- 
propriately created without any planning errors, including 
wrong source position, wrong decay correction of source 
strength and inappropriate optimization or use of an un- 
planned size or combination of applicators. However, 24 of 
169 cases had a >5% difference in total dwell time. To 
explain these differences, we investigated the correlation 
between differences in total dwell times and dh, dt and dw 
between the benchmark and individual plans. Figure 5a 
shows the relationship of differences in total dwell time 



(vertical axis) with those in dh (horizontal axis). Spearman's 
rank correlation coefficient (r s = 0.836, P <0.01) showed a 
strong relationship between the discrepancy in total dwell 
time and those in dh. In contrast, no correlations were found 
between the discrepancy in total dwell time and those in dt 
(/•„ = 0.371, P =0.075) or dw (/•„ = 0.290, P =0.149) 
(Fig. 5b and c). 

Tandem-cylinder cases 

Figure 6a shows differences in total dwell time between 
the benchmark and individual treatment plans. Differences 
averaged -1.5% ±3.1% (range, -13.0% to +0.4%), with 2 
of 29 cases exceeding 11%. 

Figure 6b shows the differences in dh fronl , dh ext , dt and 
dw between the benchmark and individual treatment plans. 
Differences averaged -1.5 mm ±4.9 mm (range, -19.0 mm 
to +4.0 mm), +1.8 mm ±5.2 mm (range, -2.2 mm to 
+20.8 mm), +0.1 mm±1.0mm (range, -1.3 mm 
to +4.3 mm) and +0.2 mm ±0.8 mm (range, -0.4 mm to 
+3.9 mm), respectively. The differences in 2 of 29 cases, 
which also exceeded an 1 1 % difference in total dwell time, 
exceeded -19 mm and +20 mm for dh front and dh exb re- 
spectively. These cases were found to have been treated 
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Fig. 5. Correlations between the differences in total dwell time and in (a) dh, (b) dt and (c) dw. Vertical axis 
shows the % differences in total dwell time between the benchmark and individual treatment plans, and the 
horizontal axis shows the differences in dh, dt or dw in mm between the benchmark and individual treatment 
plans. 



with an unplanned tandem length, resulting in an incorrect 
setting for the source dwell positions in treatment planning. 

Determination of tolerance limit 

For tandem-ovoid cases, the tolerance level of total dwell 
time, dh, dt, and dw were -7.5% to +7.2, -6.0 mm to 
+4.1 mm, -1.8 mm to +3.4 mm and -3.0 mm to +2.8 mm, 
respectively (Fig. 7a). 

For tandem-cylinder cases, two cases were excluded 
from the determination of tolerance limits because they 
were human error-related. Tolerance limits for total dwell 
time, dh front , dh ext , dt, and dw were -2.5% to +1.1%, - 
2.6 mm to +2.3 mm, -2.2 mm to +3.0 mm, -0.9 mm to 
+0.6 mm and -0.4 mm to +0.4 mm, respectively (Fig. 7b). 

DISCUSSION 

We used benchmark plans to establish a highly simple, 
easy, and fast independent treatment planning verification 
method for high-dose-rate intracavitary brachytherapy for 
carcinoma of the uterine cervix that requires no special 
skills such as developing TG43-based in-house software. 



Despite its great simplicity, analysis of a large number of 
clinical cases showed that our method was able to detect 
human error-related planning mistakes, and to evaluate the 
quality and consistency of applicator geometry. 

The Manchester method, which was first suggested by 
Tod and Meredith in 1938, has been the most broadly used 
for the treatment of carcinomas of the uterine cervix, with 
some modifications from the original [20, 21]. They 
demonstrated an 'ideal' system in which loading patterns of 
milligrams of radium were determined based on the size of 
the tandem and ovoids to achieve as constant a dose rate at 
point A as possible, no matter what combination of applica- 
tors was used, and to ensure a suitable ratio between the 
intra-uterine and vaginal contribution [21]. This rule has 
been applied to high-dose-rate brachytherapy with an 
Ir-192 stepping source by modulating the weight of dwell 
times. In our institution, manual optimization in treatment 
planning is also based on the Manchester method with 
some modifications [19]. In addition, the applicator is set 
such that its geometry is consistent with the 'ideal' geom- 
etry specified in the benchmark plan. If there were any 
planning errors and inappropriate applicator setting, total 
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dwell time and dose distribution shapes, including dh, dt 
and dw, in individual treatment plans should agree with 
those of the benchmark plan. From these points, we estab- 
lished a method for the independent verification of patient- 
specific treatment planning QA by comparing benchmark 
plans with individual treatment plans. 

Several other independent verification methods for indi- 
vidual treatment planning have been reported. Kumar et al. 
developed an in-house application that calculates the dose 
at arbitral points [13]. Lachaine et al. also developed an 
in-house application that achieves very fast calculation of 
point dose [14]. Such kind of applications are likely to be 
useful in the commissioning of treatment planning systems 
and partly also in individual treatment planning QA in 
terms of parameters such as source strength, treatment 
date and source table, which users input by themselves. 
However, because these applications use the same Cartesian 
coordinate system, digitized applicator paths and dose point 
coordinates as those in the treatment planning system, they 
are unable to detect human errors associated with the 



0) s 




Total dwell dh dt dw 
time(%) (mm) (mm) (mm) 



treatment planning process, such as setting of prescription 
point (Point A) with the wrong coordinate, the incorrect 
digitization of applicators, incorrect dose points or applica- 
tor points, improper magnification of simulation films, or 
use of an unplanned size or combination of applicators. 

Several groups have previously proposed a method of 
checking total dwell time as a fraction of treatment length 
and dose prescription, or dose area index for planar 
implants and single catheter interstitial brachytherapy 
[8-11]. Recently, for example, Das et al. reported that total 
dweU time can be predicted from the reference volume 
covered with the prescribed dose (V100) in both single cath- 
eter and multiple catheter interstitial implants by the retro- 
spective analysis of VI 00 from many clinical cases [12]. All 
these reports were focused on interstitial implants rather than 
intracavitary brachytherapy, however, and little information 
on intracavitary brachytherapy for carcinoma of the uterine 
cervix is available. In 1992, Thomadsen et al. demonstrated 
a method that assures the consistency of dwell times and 
dose distribution with previous treatment fractions and 




Total dwell dhl'ronl dhcxt dt dw 
timc(%) (mm) (mm) (mm) (mm) 



Fig. 7. Tolerance levels of (a) tandem-ovoid and (b) tandem-cylinder. 
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previous patients [7]. Although our method is basically 
similar to their concepts, we created benchmark plans in 
which the 'ideal' geometry of the tandem-ovoid or tandem- 
cylinder can be achieved because these applications are con- 
structed outside the patient's body. Therefore the geometry 
of the applicator can be evaluated in every treatment by 
comparison with that of the 'ideal' geometry in the bench- 
mark plan. 

Although many independent verification methods have 
been reported, as described above, our present study is one 
of only a few to evaluate the usefulness of the method in a 
large number of clinical cases of intracavitary brachyther- 
apy. In tandem-cylinder cases, two cases were found to 
have >11% differences in total dwell time between the 
benchmark and individual treatment plans, and >18 mm 
differences in dh &ont and dh ext . Review of these two cases 
showed that these differences were due to the unplanned 
use of tandem length, which resulted in the use of incorrect 
settings for the source dwell positions in treatment plan- 
ning. The results clearly demonstrate that our method can 
easily identify such kinds of human error. 

Regarding tandem-ovoid cases, a thorough review 
revealed no human-related planning errors. We next exam- 
ined the reason why 24 cases of tandem-ovoid cases had 
>5% differences in total dwell time between the benchmark 
and individual plans. We found that these differences 
strongly correlated with differences in dh (Fig. 6), indicat- 
ing that when the ovoid position shifts cranially compared 
with the benchmark plan, total dwell time decreases 
because the distance between the sources in the ovoids and 
point A becomes shorter. Conversely, if the ovoid position 
shifts caudally to the tandem flange compared with the 
benchmark plan, total dwell time increases, because the dis- 
tance between sources in the ovoids and point A becomes 
larger. These facts indicate that our method is useful in not 
only finding human errors and software bugs but also in 
evaluating the quality of the applicator insertion technique. 
In other words, the evaluation of both the differences in 
total dwell time and dh could provide a good indicator for 
the quality of the applicator's geometry. 

We set tolerance limits for differences in total dwell time 
and dose shape between the benchmark and individual 
treatment plans. Ezell et al. reported that they set action 
levels for per-patient intensity modulated radiation therapy 
verification using confidence limits [25], as first proposed 
by Venselaar et al. [24]. If the confidence limit is estab- 
lished with sufficient points to provide good statistics, then 
the value of 1.96o suggests that variations in excesses of 
the limit would occur about 5% of the time. We determined 
the tolerance limits by using confidence limits for total 
dwell time, dt, dh (dh front , and dh ext for tandem-cylinder 
cases) and dw. To calculate tolerance limits, we excluded 
the two cases with >18% differences in tandem-cylinder 
cases to eliminate the effect of human error-related planning 



mistakes. These tolerance limits might be one indicator in 
the evaluation of individual treatment plans (i.e. rechecking 
of treatment planning, use of appropriate size of applicators, 
inappropriate applicator geometry, etc.). 

One limitation of our study warrants mention, namely 
that our method is useful only for Manchester-based intra- 
cavitary brachytherapy. For carcinoma of the uterus, 
however, most treatment centers in the world have followed 
a traditional concept based on the Manchester loading pat- 
terns [3, 26]. Moreover, local control rates are significantly 
improved with EBRT and brachytherapy based on the 
Manchester method [1]. Our method therefore appears 
useful despite this limitation. 

In conclusion, we established a highly simple, easy and 
quick independent verification method using benchmark 
plans for intracavitary brachytherapy based on the 
2D-based Manchester method. Despite the great simplicity, 
our method can evaluate the quality of the applicator inser- 
tion technique, as well as identify human errors in treat- 
ment planning. 
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